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Abstract 

Consider n non-intersecting Brownian particles on K (Dyson Brownian mo- 
tions), all starting from the origin at time f = 0, and forced to return to jc = 
at time f = 1 . For large n, the average mean density of particles has its support, 
for each < f < 1, on the interval ±\/2n/(l — t). The Airy process ^{t) is de- 
fined as the motion of these non-intersecting Brownian motions for large n, but 
viewed from the curve ^ : y = sj2nt{\ — t) with an appropriate space-time 
rescaling. Assume now a finite number r of these particles are forced to a differ- 
ent target point, say a = Po\/«/2 > 0. Does it affect the Brownian fluctuations 
along the curve 'ta for large nl In this paper, we show that no new process ap- 
pears as long as one considers points {y,t) e 'lo, such that 0</<(l+Pg)^', 
which is the /-coordinate of the point of tangency of the tangent to the curve 
passing through {pQ^/nJl, 1). At this point of tangency the fluctuations obey a 
new statistics, which we call the Airy process with r outliers s^^''^ (t) (in short: 
r-Airy process). The log of the probability that at time x none of the particles in 
the cloud exceeds x is given by the Fredholm determinant of a new kernel (ex- 
tending the Airy kernel) and it satisfies a non-linear PDF in x and T, from which 
the asymptotic behavior of the process can be deduced for T — oo. This kernel 
is closely related to one found by Baik, Ben Arous and Peche in the context of 
multivariate statistics. 
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Dyson |[T6l made the important observation that putting dynamics into random 
matrix models leads to finitely many non-intersecting Brownian motions (on M) for 
the eigenvalues. Applying scaling limits to the random matrix models, combined 
with Dyson's dynamics, then leads to interesting infinitely many diffusions for the 
eigenvalues. This paper studies a model, which stems from multivariate statistics 
and which interpolates between the Airy and Pearcey processes. 

Consider n non-intersecting Brownian particles on the real line M, 

-oo<Xi(0 <-..<Xn{t) <oo, 

with (local) Brownian transition probability given by 

1 (v--v)2 

(0.1) p{t;x,y) := —=e < , 

all starting from the origin x = at time t = 0, and forced to return to = at time 
t = I. For very large n, the average mean density of particles has its support, for 
each < f < 1, on the interval {—y^2nt{l —t), yj2nt[\ —t)), as sketched in Figure 
[OB 

The Airy process A(t) is defined as the motion of these non-intersecting Brow- 
nian motions for large n, but viewed from an observer on the (right hand) edge- 
curve 

^ : {y= yj2nt{l-t) > 0} 

of the set of particles, with space stretched by the customary GUE-edge rescaling 
n^/^ and time rescaled by the factor n'/^ in tune with the Brownian motion space- 
time rescaling; this is to say that in this new scale, slowed down microscopically, 
the left-most particles appear infinitely far and the time horizon t = 1 lies in the very 
remote future. Thus, the Airy process describes the fluctuations of the Brownian 
particles near the edge-curve looked at through a magnifying glass, as shown in 



Figure 0. 1 



The Airy process was introduced by Prahofer and Spohn ESI and further inves- 
tigated in |[20l [211 [27] [3. Notice that in this work the Airy process is not viewed 
as the motion of the largest particle (point process), but as the motion of the cloud 
of particles, which will be described as a determinantal process. Giving a pathwise 
description of this motion remains an open problem. 
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Figure 0. 1 . Airy process 



Assume now that, among those n paths, < r < « are forced to reach a given 
final target a > 0, wliile the {n — r) remaining particles return to the position x = 0. 
Consider the probability that, at any given time < f < 1 , all of the particles avoid 
a window £" C M, namely (the superscript in the probability 
refers to the target points) 
(0.2) 

all;cy(0) = 



here and later 



jOfl 
Br 



{all Xj{t) eE' 



all Xj{t) G 



r right paths end up at a at f = 1 
{n — r) paths end up at at ? = 1 



Does the fact that a finite number r of particles are forced to a different target 
point, in particular the target point a = po\/i/2 > for some arbitrary parameter 
Po, affect the Brownian fluctuations along the curve for very large n? It is 
understood here that near the points of the curve (under consideration), one uses 
the same scaling as the Airy process. In this paper, we show that no new process 
appears as long as one considers points 



(0.3) 



{y,t) e such that 0<t < 



1 



1+Po 
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Figure 0.2. The r-Airy process 



Observe that t = to = { \ + p^) ^ corresponds to the point of tangency (^0,^0) of 
the tangent to the curve passing through the point {po\/n/2, 1) of the {t = l)-axis; 
i.e. 



(0.4) {yo,to) = 1 1 e ^. 



P()V2w l_ 

l+Po''l+Po 



At this point of tangency (yo,?()) the fluctuations obey a new statistics, which we 
call the Airy process with r outliers (t) (in short: r-Airy process); see Figure 
0.2 In particular, for the target point a = ^Jnjl, the r-Airy process occurs at the 
the maximum of the edge curve i.e., po = 1 and thus fo = 1/2. Notice, the r- 
Airy process is an extension of the Airy process; the r-Airy process coincides with 
the Airy process, when the target point a coincides with or, what is the same, 
when r = 0. The Airy process is stationary, whereas the r-Airy process ceases to 
be stationary. 

Given a target point a = po \/n/2, the point of tangency of the tangent to the 



curve ^, passing through a, can be written, in accordance with (0.4 1, 



(0.5) {yQ,to) = — — 2 > — 2 = (3^0, ?o = ^ , ^ e 'T, 

\ 1+Po 1+Po / \ 2coshTo 2coshTo / 
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upon introducing a new parameter Tq, whose significance as a new time parameter 
will be clear from section [T] (formula ( |1.8| )), and wliicli is defined by 



(0.6) e-'°:=po 



Given a subset £ C M, the r-Airy process in the new time T will be obtained by 
taking the following limit: 



(0.7) P(j^W(T)n£ = 0) 



2n + - 



V ' Vl+^"'^'"^"^' J 2cosh(To+;^) ; ■ 

Notice that for T = and upon ignoring the set E'^, space and time on the right hand 
side of the formula above equal {yo,to) as in (0.5 1. In section|2j the limit above will 



be shown to exist, independently of the parameter po, at the same time establishing 
some universality. 

1 0J_ below, define the function^ 



In order to state Theorem 



(0.8) Af{u; t) := / {^ia - ^ 

Jc 27t 



and the standard Airy function A(m) := Aq = AJ" satisfying the differ- 
ential equation A"(x) = xA(x). Given these functions, the Airy and r-Airy kernels 
are defined by (|0.9[) below and will be used in Theorem |0.1[ 



(0.9) 



poo 

/^(o)(m,v)= / dwA{w + u)A{w + v) 
Jo 

k[''\u,v) = i (iwA7(w + M;T)A+(w + v;T). 

JO 



A{u)A' {v) - A' {u)A{v) 



u — v 



We now state: 

Theorem 0.1. Consider non-intersecting Brownian motions as above with r par- 
ticles forced to a target point Po\/«/2 > at time t = \. For large n, the average 
mean density of particles as a function oft has its support on a region bounded to 
the right by the curve ^to (defined just below (|0.i|)). The tangent line to passing 



through {po^yn/2,l) has its point of tangency at {yo,to), given by (0.41. Letting 
n — > oo and given an arbitrary point {y,t) G ^ with t <tQ, a phase transition occurs 
at t = tQ, which is conveniently expressed in terms of the new parameters defined 



' C is a contour running from oog^'^r/fi jq ^gin/b^ f,\xch that — /T lies above tiie contour. 
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Figure 0.3. Pearcey process 



b^t = {\+e ' and to = {l+e ^^") ^; one has the following limit: 



(0.10) Yim¥t''^""-Hallxi 



2n + 



2cosh(a + -^; 



T)n£ = 0) =det(/-i«:(''))^ for 0<t<tQ 

( i.e., < a < To), 

'•)(T)n£ = 0) =det(^/-4''^)^ for t = to 

( i.e., a = To). 

These Fredholm determinants are genuine probability distributions for each value 
of Z. Notice that for < a < To, the limit above is independent of the time T and 
the number of outliers r, unlike the case o = Tq. 



As is clear from Figure 0.2 and for t > tg, the r outliers separate from the bulk 
as a group; thus after an appropriate shift to take into account their displacement 
from the curve and using a different scaling, they behave probabilistically like 
the r eigenvalues of an r x r matrix in the GUE ensemble. 

The r-Airy process can also be viewed as an interpolation between the Airy 
and Pearcey processes. The Pearcey process Il28ll23l l4ll5l is defined as the hmit of 



Remember from 1 0.6 1, one has 



Po ■■ 
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non-intersecting Brownian motions (for large n), all leaving from at time t = 0, 
with {I — p)n paths forced to end up at and pn paths forced to end up at 
at time t = I; see Figure 0.3. The Pearcey process then describes this Brownian 
cloud of particles (for n ^ oo) near the time t of bifurcation, where the support of 
the average mean density goes from one interval into two intervals, with stretched 
space and time. Then the boundary of the support of the average mean density of 
particles in {y,t)-spsice has a cusp. In sectionjsjit will be shown sketchily how the 
location of the cusp, when the proportion p of particles tends to as r/n, tends to 
the precise place {yo,to) where the r-Airy process occurs. 

Remark 0.2. To simplify the notation, one often writes 



Remark 0.3. The joint probabilities for the r-Airy processes for a finite number of 
times can be defined in a similar way and lead to a matrix Fredholm determinant, 
to be discussed in a later paper. 

This phenomenon is closely related to statistical work by Baik-Ben Arous- 
Peche im, Baik ||9l and Peche ll29l . Indeed, consider a (complex) Gaussian popu- 
lation ~y G C^, with covariance matrix £. Given M samples yi, . . . ,yM, the (cen- 
tered) sample covariance matrix S := -hXX^ , where 



is a positive definite matrix and is an estimator of the true covariance matrix £. 
One may test the statistical hypothesis that £ = / or that £ has all eigenvalues 
= I, except for a few outliers. When all the eigenvalues of £ are = 1, then the 
limit distribution of the largest eigenvalue of the sample covariance matrix S, for 
N and M tending to oo in the same way, is given by the Tracy- Widom distribution. 
Then Baik-BenArous-Peche [8 ] noticed that, this is still so, if the eigenvalues of £ 
contain some outliers, which are not too large. There is a critical point at which 
and beyond which the largest eigenvalue of the sample covariance matrix S will 
be different from the Tracy- Widom distribution. At this point of phase transition 
Baik-BenArous-Peche JH have found an Airy-type distribution in x, which is given 
by the Fredholm determinant 



where r denotes the number of eigenvalues of £ that are equal to 1 + 7 , while 
all the others are = 1; 7 is such that M/N = for M and N very large. This 

distribution was further generalized in |[8l to the case where T 7^ in the kernel k'^''^ . 
Baik proved in [9] that the Fredholm determinant of K^^^ is a genuine probability 
distribution. In the statistical problem above, the covariance matrix S is positive 
definite and therefore its eigenvalues satisfy Laguerre-type distributions; this idea 



(0.11) 





det(/-4''0 
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was extended to GUE-type distributions by S. Peche (29]. In the present paper one 
finds that the shift T appearing in the kernel K^:' ' is precisely the rescaled time of the 
non-intersecting Brownian motion model! The arguments will appear in sections [T] 
and 12 

This paper also shows that the probability of the r-Airy process or, what is the 
same, the Fredholm determinant det ( / — Kif"^ ) satisfies a non-linear PDE in x 



(jr,~) 

and T, depending on the number r of outliers, as established in section[6| 
Theorem 0.4. The logarithm of the probability 

Q{x,x) :=logP(supj#W(T) <x) =logdet(/-i«:i''^)^^^^ 

satisfies the followins non-linear PD^^ with both, the fimction 2(t,x) and the 
PDE, being invarianvjunder the involution (t,x, r) (— T,x, — r), 




dtdx J \dx^ dtdx J dl dxdl ^ dt^ 



- '3^Q ^yQ_3^' 



dz^ dzdx dx^ 

with "initial condition" , given by the log of the Tracy-Widom distribution, 

poo 

(0.13) Qo{x) := lim G(t,x) := logP(supi^(T) <x) = - {a - x)g^{a)da, 

Jx 

where g{(x) is the Hastings-MacLeod solution ofPainleve II, 

(0.14) g" = ag + 2g\ withg{a)^ for a /oo. 

Remark 0.5. Obviously, the PDE ( 0.12[ ) has the following structure 
(0 15) ^(^\^d^_^lQd^ 

Remark 0.6. The following simple recipe gives the PDE f or ¥{.s/^''\t) = 0) 
for a general set £ = Uj [x2,-i,X2,], replacing the PDE (0. 12i for £" = (x,oo); indeed, 
perform the replacements 

/t-i 




^ The Wronskian {f,g]x with regard to the variable x is defined as f'g — fg' 
^The in 
Lemma 7.5 



^The invariance under the involution is obvious for the equation; for the function Q{'C,x), see 
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with the understanding that differentiation must always be pulled to the right. 



Although the average mean density of the particles is insensitive to the pres- 
ence of outliers, the presence of the r particles forced to reach the target a > at 
f = 1 is already felt, when ? ^ in the f-scale; that is when T — > — oo in the T- 
scale. The net effect is that it pulls the edge of the cloud of particles in the average 
towards the right to first order like |r/T|; i.e., the more so when r gets large; this 
edge then behaves like the Airy process shifted in space, up to and including order 
1 / T^.The PDE in Theorem 0.4 is a convenient instrument to extract the remote past 



asymptotics, as shown in Section[7]and stated in the theorem below. 



Theorem 0.7. The PDE with the initial condition Qo{x) as in (0.131 admits the 
asymptotic solution, for T —oo, of the form 



Qo 



r s. , r ^ r 



For the probability itself, one haQ 
P(sup^('''(T) <;c) 



sup^(T)<(x+^)(l + ^ 



1 + 



5T' 



For T 



-oo, the mean and variance of the right edge of the process behave as 



E(sup^W(T)) = E(sup^(°)(T)) ( 1 



3t3 



var (sup^W(T)) = var(supi^(")(T)) I 1 



lr\ ,\ 



dMe-^_-YAM^-2AM) 



The probability ( |0.2[ ) is related, via a change of variables, to a Gaussian matrix 
model with external potential, 

-/ 

where M'n(E) denotes the Hermitian matrices with eigenvalues in £ C M and with 
a diagonal matrix 

/ a 



(0.16) 



a 



O 



O 



/ 



with ^5 :=x'e;, +4x20 + + lo/r eo -6j;'dyrduQ\ 



,1/2 
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The relationship between the Karlin-McGregor non-intersecting Brownian motions 
and matrix models has been developped by Johansson [19] and, for the Gaussian 
matrix model with external potential, by Aptekarev-Bleher-Kuijlaars [7 |. The latter 
model has come up in many other situations; among others, see ll24l[miT2l|T3ll32l 
[3l|20l Uni im in. Following the method of Adler-van Moerbeke f4\, using the 
multi-component KP hierarchy (section[3]l and the Virasoro constraints (section[4]), 
it is shown in section [5] that this matrix model satisfies a non-linear PDE in a and 
the boundary points of E. The PDE of Theorem 0.4 then follows from making an 
asymptotic analysis on that PDE. The asymptotic behavior of the r-Airy process 
for T —oo follows from solving the partial differential equation near T = —oo 
for the Airy-kernel initial condition; this process introduces some free constants, 
which then can be determined by the asymptotic properties of the r-Airy kernel for 
X ^ oo. It is an interesting question whether these results can be deduced via the 
Riemann-Hilbert methods, in the style of JTSll . 



1 A constrained Brownian motion with a few outliers 



The purpose of this section is to motivate the scaling limit (0.7 1 leading to 
the definition of the r-Airy process. The Airy process was originally defined as 
an edge scaling limit of Dyson's non-intersecting Brownian motions, in the same 
way that the Tracy-Widom distribution was obtained as an edge scaling limit of 



QUE. The non-intersecting Brownian motions (0.2 1 with target point a = can be 
transformed to the Dyson Brownian motion, as will be explained in this section. 
This transformation will be used in the case of a a 7^ target. In the course of doing 
this, one must consider a Gaussian matrix model with external potential. 



1.1 Dyson Brownian motion 

The Dyson process ( 11161 ) describes the motion of the eigenvalues of an 
nxn Hermitian random matrix B{t') whose real and imaginary entries perform in- 
dependent Ornstein-Uhlenbeck processes, given that the initial distribution is given 
by invariant measure for the process, namely 

Then the process is stationary and its probability distribution at any time is given 
by 

(1.1) ¥D^.{^\\Xi{t')eE')=Z-^ I dBe-^'^' = Z-^ [ dMe-^-'^'"'. 

Jje„(E') Jj^„{V2E') 
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The probability for two times < t[ < t2 is then given by (set c' = e 



P(all (B(f^) -eigenvalues) G E[, all (B(f2)-eigenvalues) G E2) 



(1.2) 



all B, -eigenvalues eE[ / ^ _ ^l2\n- /2 

all B2 -eigenvalues ^ ' 



Z'-'dM,dM2 e-iTr(MKM|-2c'M,M,) 



1.2 Constrained Brownian motion with target a > 

As in the introduction, consider the n non-intersecting Brownian particles on 
M, all starting from the origin at time t = 0, where among those paths, 1 < r < n 
are forced to end up at the target a > 0, while the {n — r) other paths return to the 
position ;c = at time t = 1. Remember Fg"^ denotes the transition probability. Then 
the probability that all the particles belong to some window £ C M at a given time 
< f < 1, can be expressed in terms of a Gaussian matrix model with an external 
potential, using the Karlin-McGregor formula Il22ll for non-intersecting Brownian 
motions (see ifTTl IT9l [7] IH), itself involving the Brownian transition probability 
(|0.1|), namely: 



all Xj{t) G E 



lim 



all Yi^O JE"Zn{y,5) 1 

5i , . . . , 5n-r 
5,1-r+l ) • • • ) 5„ — > a 



all;cy(0) = 
r right paths end up at a at ? = 1 
{n — r) paths end up at at f = 1 

1 " 

Y\dxi 



det(7?(?;7;-,xy))j<. .<„det(/7(l-f;;c,v,5/))i<,vj/<„, 
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1 
1 



^M^'-jTr(M2-2AM) 



(=1 



(=1 



with A as in ( |0.16| ). In sliort, the conditioned Brownian motion is related to a 
Gaussian matrix model with an external potential A as follows: 



(1.3) ¥%{&\\xj{t)(^E)=¥„{a,E) with E = E. — 



(1-0 



, a = a 



2t 



The joint probability for the constrained Brownian motion at two times is related 
to a chain of two Gaussian matrix models (see 121) with an external potential A, 
which again by Karlin-McGregor reads as follows: 

(1.4) ¥l'l{sillxi{h)eEuallxi{t2)eE2)=K{a;c';Ei,E2), 

where (A is the same diagonal matrix as in ( |0.16| )) 

Zn Jj&,(E,)xX(Eo) 



with 



Ei=Ei 



2t2 



{t2-tl)h 



, E2—E2\ 



{l-t2){t2-h) 



(1.5) 



'\-t2)ti Vl-'l 



;i-fih 



a = a\ 



2{t2-h] 



'\-t2){\-h) 



1.3 Comparing the Dyson and constrained Brownian motions with tar- 
get a = 

From the identities ( |1.1[ ) and ( 1.3 1, one deduces the following set of identities 
for one time t' (in the time t' of the Dyson process), setting the point a = 0, 



(1.6) 
where 



Vv(all Xj{t')eE')=¥n{a,E) 
y/2E' = E = E. 



-t{^\\xj{t)&E) 



a=0 



a=0 



t{\-t) 



A,.(x) denotes the Vandermonde determinant, with regard to the r variables x = (xj , . . . ,Xr) and 
A„(x,>') denotes the Vandermonde determinant, with regard to the n variables x = (xi , . . . ,Xr) and 

y = {yi.---.yn-r)- 
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Similarly, for two times t[ and t!^, one deduces from the identities (1.2 1 and ( 1.4i, 
PDy(all Xi{t[) e E[, all Xi{t!^) e E'^) = P„ (a;c';£i,£2) 



a=0 



¥X{&\\Xi{ti)eEi,&\\Xi{t2)eE2) 



a=0 



where, as follows from ( |1.2[ ) and one deduces 



2t2 



(1.7) 
with 



(?2-fl)fl 



E2=E2] 



2{\-h) 



[l-t2){t2-h] 



\-ti 



J2_ 
l-f2 



Thus the processes are related by a clock change 

1 



(1.8) 



ti(\-ti)- 



ti 



l+e 



-2t 



2coshf-, 



and thus 1 — c'^ = . '[ . . Comparing extremities in ( 1.7 1, one finds 



E[ 



2t2 



{t2-h)ti 



l-c'2)/2 = £i 



1 



E' 



V (l-?2)(?2-?l) 

This fact, combined with ( |1.8| ), yields 
(1.9) 



■\-c'^)/2 = E2 



^/h{^\--h) 

1 



^t2{\-t2)' 



Ei = E[^ti{\-ti) 



e: 



2 cosh 



So, summarizing ( 1.8 1 and ( 1.9 1, one has the relation between the parameters of 
the Dyson process and constrained Brownian motions. 



(1.10) 



ti 



ti{l-ti) 



V~2E[. 



Theorem 1.1. (Tracy-Widom II26II , Adler-van Moerbeke ||3lj Taking a limit on the 
Dyson process, in an appropriate time and space scale, one finds the Airy process, 
which is stationary: 



r)<x)=F{x), 



Urn Foy all A, 



zV3 



G {—oo^^/2n + 



V2«V6- 



14 



M. ADLER, J. DELEPINE, AND P. VAN MOERBEKE 



F{x) being the Tracy-Widom distribution. Similarly the limit of the joint probability 
for the Dyson process yields the joint probability for the Airy process: 



lim Fdv 



2n + 



Its logarithm (setting s= (T2 — Ti)/2j 

H[s-x,y) :=logP(i/(Ti) <x + y,. 
satisfies the Airy PDF 



{T:2)<x-y), 



(1.11) 



Is 



dsdxdy 



2s' 



dy 



dy^ 



+ 



\ dxdy ' dx^ 



Corollary 1.2. Taking an appropriate scaling limit on the constrained Brownian 
motion, we have 



lim 



Br 



all xi 



1 



2n + 



a=0 



T)<X) 



l+g-2T/«i/3y 2cosh(T/«i/3) 
and similarly for two times. 

Proof. This statement follows immediately from Theorem 1 . 1 and the correspon- 



dence ( |1.6| ), upon using the clock change ( |1.8| ) and the space change ( |1.9| ), with the 
appropriate time and space scalings of Theorem |l.l[ □ 



Remembering the definition 

P„(a,£) 



1 



for the diagonal matrix A as in (0. 16 1, we now state a Theorem of Peche, which is 
closely related to the multivariate statistical problem mentioned in the introduction; 
seelH. 



Theorem 1.3. (Peche f29)) 



(1.12) 



limP„ (pV^,(-oo,2V^ + 



.1/6- 



F{x)for p < 1 
FW(x)/orp = 1 



T=0 



where 

(1-13) FW(x) = det(/-4''^;^(.H(3')) 

is the Fredholm determinant of the kernel (see ( |0.8p j 

(1.14) K'(\x,y)= I duA;(x + u;x)A+{y + u;i). 

Jo 



For r = 0, (1.131 yields the Tracy-Widom distribution F{x). 
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The Airy process will now be deformed in a way which is compatible with 
Peche's Theorem, especially concerning the target point a, 



Br 



l+^,-2T/«i/3y 2cosh(T/?iV3) 



Since, from ( 1.3 1, 



-t{^\\xj{t)eE)=¥AaJ^^;E^ 



t{\-t) 



holds, we have 



Br I l^^^_2T/«l/3y - 2cOSh(T/«l/3) 

= P„ (aV2^'^/«'^';(-oo,2V« + ^; 

Comparing this formula at time f = 1/2, or what is the same at T = 0, with formula 
( 1.12 1 suggests the choice 

(1.15) " = ^\/f- 

We thus define the Airy process with r outliers (in short: r-Airy process) £!/^''\x) 
by means of the Airy scaling as in Corollary [L2] and with the choice of a above for 
p = 1. Notice from ( 1.3 1 and using ( |1.10[ ), this can also be expressed in terms of 
the matrix model with external potential: 



'-) (T) n £ = 0) := lim all x,- / ^ ^ ^^""^"'^^ 



2n+ ■ 



l+e-2T/ni/3y 2cosh(T/«V3) 



In a similar way, one defines the joint probability of s^''''\t) for any number of 
times. The next section deals with this limit expressed in terms of a Fredholm 
determinant. 



2 The existence of the Umit to the r-Airy kernel 



The first part of this section deals with a sketch of the proof of Theorem 0. 1 



details and rigor can be found in [8] and [|291 . In the second part, the r-Airy kernel 
(|1.14|) will be expanded for large time T. 



Proof of Theorem 0.1 From the explicit Brownian motion transition probability, 
one shows (see Johansson llT9l and Tracy-Widom ||28]| ) 
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'^^;(allxXOG£)=det(/-//W)^,(^.) 



where 

Hl[\x,y)dy: 
(2.1) 



dy 



27r2(i 



SI 



dz I dw e i-' + i-f + i-r — 



w — a 
z-a 



w - 



^ is a closed contour containing the points and a, which is to the left of the line 
Fl := L + /M by picking L large enough. So, 9?(w — z) > 0. Consider now an arbi- 
trary point {y,t) on the curve parametrized by (2.2 1 and the point (p a/«/2, 1), 
which is the point of intersection of the tangent to 'tf at {y,t) with the axis (f = 1); 
as pointed out before, it is convenient to parametrize p by p = and thus 



(2.2) (3^,0 



pV2?i 1 

T+p2'T+^ 



2n 



2cosh(a)' 1+£>-2ct 



Consider the Brownian motions with r outliers forced to a point pQ-^fnfl with 
< Po < p at time t = \\ here also parametrize po by po = and set a = 
Po/P 



The main issue is to compute the following limit, for < po < p 



lim 



j(0,poV"/2) 
Br 



all X; 



1 



2?i + 



l+^' 



„i/3 ' 



2cosh(a + 




which dictates the space- and time-scale to be used in the kernel for large n. 
Let t and y be the time and space variables for the Brownian motion, which in terms 
of the new time and space scale 5£ reads: 
(2.3) 



1 



In- 



2n + 



l+e-2{f7+T/«l/3)' 



2cosh(a + 



2cosh(a + ;^)' 

with target point a = po\/n/2 = e^''^°y^n/2. Putting this rescaling in the integral 
(2.1 1 suggests changes of integration variables 



z:=Z\l-e''^ ^ andw:=ww-e 



in the integral (2. 1 1; the exponential will contain a function F{z), with Taylor series 
atz = 1: 

(2.4) = ^ - 2z + logz = F(l) + -(z- 1)3 + 1)1 



Also set ^ := 1 + for some parameter 7 and 



(2.5) 



(l-a)'>"'^(') for 0<a<l 
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Then, using, 



one check^ using the rescaling (2.3 1 



-yn'/^ (w-z) 



dy, for 9^(^-1) > 0, 



H!^''\x,y)dy\ 



2% I Jsi Jtl 



w \ I w — ae 
z J \ w 



z-oce 



1 



J «l/3(w)-z) 

271 Jcj Vz-ae^/"^ 



-e 

2nZ„ jyl 



yl\u+y){z-i;) 



- f dw 

■n JTl 



W 



w — ae 



t/«1/3 



,-«l/3(,+y)(,p_Q 



Conjugating the kernel, which leaves invariant the Fredholm determinant, one finds 



«l/3(,_„) (r) 



H!f^ {x,y)\^dy = -dve^'^"-'-^ l^^dy^i"\u+y)^i"\v + y) 



with 



(2.6) 



I dze-''^^^ 
J9 



271 Zn JYl 



z-ae 



w 



t/«1/3 



w — ae 



t/«1/3 



by elementary computation 



- 7 

t 2(cT+-rW) dy dzdw n^dzdwdv 2xz ,^ i,. tz n __2 

-, =e , = ; — :— , "i = {2n + un^)z, = • 

l-t 1-fvv-z w-z 1-f 1-? 2 
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Using the Taylor series (2.4 1 for f (z) and the value ( 2.5 1 for Z„, one is led naturally 
to pick a new variable u such that z — 1 



nl/3 



. Then 



\z-ae^/" ' J 

.„./3,,-(.-„.„-*,o(-i,)) 



u + n 



1/3 



M + «l/3_ah+-^ + .. Ul/3 



,.3 



e ^'^du e "3 (;^)'^ + lower order terms for a = 1 
e^'^^Jw + lower order terms for 0<a<l. 



Upon remembering the definition of the functions Af{x,z) and the Airy function 
A(x) and also the fact that a = po/p, the rigorous saddle point argument, given in 
IBi^, yields 



lim^i"'(x) 



(2.7) 



(2.8) 



271 



J_ 

271 



for po = p 
for < Po < p 



lim A'"'(3') 



for Po = p 
for < Po < p 



where 'S' is a contour running from ooe^'^l^ to ooe^'^l'^, with an indentation to the 
right of T , such that T is to the left of the contour, and where is a contour 
running from oae^'^l^ to ooe^^l^ , with an indentation to the right of T, such that T 
lies also to the left of the contour. Upon rotating the two contours and deforming 
slightly, since the integrand is pole-free, one gets the final identities in the 
equations above. 
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Therefore P(sup i;/^'") (r) n £ = 0) = det(/ - k'(^ )e , with 

Ki''\u,v)dv = Um e"'"^''-"^H!f^ {x,y) dy 

dv Jq dw {u + w, t)A+(v + w, t) for Po = p 

/o°° A(m + w, t)A(v + w, t) for < Po < P . 

Baik [9] has shown that the Fredholm determinant of the r-Airy kernel is a proba- 
biUty distribution, i.e., 



lim detH-Ki''^ 



1 

-r^±~ V " J (.v,oo) 1^ 0. 

This estabUshes Theorem lO.il □ 

Remark 2. 1. In this section one lets T —>■ — oo, which implies that —iz remains above 
the contour C and is thus compatible with the contour mentioned above. Letting 
z — > +00 would require a drastic change of the functions Af. 

The next statement concerns the asymptotic behavior of the r-Airy kernel for 

(2.9) Ki''\u,v)= [ dwA;{u + w;T)A+{v + w;T), 

Jo 

where (remember) 

(2.10) Af{u;T)= [ 

Jc 271 

where C is a contour running from ooe^'"^/^ to ooe"^/^, such that —ir lies above the 
contour. This limit is compatible with the contour C appearing in the definition of 
the functions Af, since then —ix remains above the contour C, as required. 

Lemma 2.2. Given the "initial condition " 

M, ^ rnw ^ Aiu)A' iv) - A' iu)Aiv) 
Km K\'Hu,v) = K^°\u,v) := ^ ^ ^ ^ y ) y I ^ "Airy kernel", 

U — V 

the kernel K^''\u,v) behaves asymptotically for T — > as 

i^{r) Ar) j^(r) 

k(^)(u,v)=Ko + ^ + \ + \ + ..., 
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wherd^ [for n> \, 



n-l 



n\ \ dw 
'2{n-2)\ \dw 



A{u + w)A{v + w) 



w=0 



n-2 



A'{u + w)A(v + w) -A{u + w)A'(v + w) 



(«-3)! 



V 



3n-l / d 



24 \ 5w 



n-l 



A(m + h')A(v + w) 



«-3 



"-i^l-l A'(« + w)A'(v + w) 



2 V5 

{polynomial of degree n — 2> in r). 



w=0 



Although the kernel K\ (m,v) involves integration, the terms K- (m,v) m the ex- 
pansion never involve integration, they are quadratic in the Airy function and its 
derivatives; also the Kj''\u,v) are polynomials in r of degree i, divisible by r, with 
alternately symmetric and skew-symmetric coefficients in u and v, the top coeffi- 
cient being symmetric. 



Proof. In order to expand the kernel (0.9 1 with regard to T for T ^ — oo, set the 
expressions ( |2.10[ ) into the kernel (0.9 1, which then becomes a triple integral. Set 
a = ia and j3 = ib and consider the following Taylor expansions about T = — oo, 

/ a+l \ 



1 



(a-T)'-(-j8 -t)-'" 



1 + 



1 



=» 1 ["2'] 

n=\ ^ j=0 



(r(a + /3))"-2ve2;(«,j8) ^ 
+(r(a + ^)r-2i-i(a-/3)(22;(a,i8), 



(2.11) 



1 



1 + 



(K«+]B))" I (r(a+P))"-'(«-;8) 
11! 2(«-2)! 
(ria+P))"-^ f {3n-l)(a+pf (n-l)ap 



+ 



(n-3)! 
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+ .. 



where Q2j and Q2j are symmetric homogeneous polynomials of degree 2j in the 
arguments, since the first expression is invariant under the involution r — r and 
a ^ The coefficients of 1/t" are divisible by r, for the simple reason that for 
r = 0, the expression above equals 1. 

Also notice multiplication by ia of the integrand in the kernel (0.9 1, 



Ki''\u,v)= I dwf e-3 
dw 



1 



la — T 



da 
In 



/,ih^+lh{w+v)^_.j^_^ydb 



In 



4n^ T 



Whenever d/dw appears with a negative exponent in the formula below, it is set = 0. 
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can be realized by taking and similarly multiplication of the integrand by ib 
is realized by taking ^^^^ ; thus we have the following recipe 



la 



ib 



du 
dv 



dwA' {u+w)A{y+w) 
dwA{u+w)A' iy+w) , 



and so in particular, 

{iaf'iibf^iia+ibf ^ 



d d 
du dv 



d 

du 



A{u+w)A{y+w)dw 



n-l 



a('=')(m + w)a('=2)(v + 



w 



w=0 



Notice that, since ia + ib factors out of every term in the expansion ( |2.1 1[ ), the 
kernels obtained never contain integration. In addition, since ia — ib factors out 
of every other term, every other term in Kf\u,v) must be skew; in particular it 

vanishes for u = v. One then reads off the K^''^ (m, v) 's from the expansion (2.111 and 
the recipe above, upon using occasionally the differential equation A"{x) = xA{x) 
for the Airy function, thus ending the proof of Lemma [2!2| □ 



Remark 2.3. As an example, we give explicit expressions for the first few K^' \u,v) 's: 
(2.12) 

—r A{u)A{v) 



Ki;\u,v) 



u — v 



k['-\u,v) 



4\u,v) 
4'\u,v) 



'-{A'{u)A{v) +A{u)A'{v)) + '-{u - v)K^'>Hu,v) 

3 2 

-'-{A"{u)A{v) + 2A'{u)A'{v) +A{u)A"{v)) + "-{v - u)A{u)A{v) 
6 2 

{A"iu)Aiv) +A{u)A"{v) -A'iu)A'{v)) . 



In order to find the PDE for the transition probability, one will need an estimate 
on how the actual transition probability for the finite problem converges for « — > oo. 
This will be used in (6.3 1. 

Corollary 2.4. Forx G M sufficiently large, one has for some constant C > 0, 



(2.13) 



1 Tn>(°'''o 



all Xi 



1 



2n + 



-2(To-+ 



< 



,1/3' 



2cosh(To+^) 



-logP(supi/W(T) <x) <Cn-^^\ 
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Proof. For any trace class operators Kn and K^o, set[^ 

M:=max(||^„||,||^„||). 
Then assuming M < 1, one checks, using K]^ —Kl, = Y,'fLQKi,{Kn — Koo)K',^^^^ and 



Blljj. and \\AB\\^^ < 



tri 



(2.14) |log det(/ -/:„)- log det(/-^ec)| = |tr(log(/-/:„) -log(/-^„))| 

oo 



ly-i ly-i 

tr( " . °° ) 



< Lyll^i-^^lltr 
1 ' 

oo 



< 



IKf, — Koc 



'tr 



Setting (see notation ( |2.6| )) 

K„{u,v) := 



dyj^i"\u+y)^r{v + y) 

poo 

K„iu,v) := e^'-"'^^ / dyA;{u + y,T)A+{v + y,T), 
Jo 

one checks 

\Kn{u,v) - Koo{u,v)\ 

< £dy(^J^i''\u+y)-e-^^'+y^U;{u+y,T))ji"\v + y) 

dye-^"+y^U;{u + y, T) (^i"^(v +3;) - ^'('■+>')M+(v + J, T) 

An argument similar to the one of Baik-BenArous-Peche lH |29l shows that for 
given xq € M, there are constants C > 0, A'^ > 0, 



l-M 







r(«) 



n 



1/3 



, for « > A/^ and x > xq. 



Define the three norms on a Hilbert space: the sup, the trace and the Hilbert-Schmidt norms, 

with||r||<||r|L<||r||j^: 

||r||= sup^ = sup|A,i ||r|ltr = Tr(rr)i/2 = j:iA,.|, = (Trrr)V2 = (j^|i,|2)i/2. 
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Viewing the functions in the integrals on the right hand side as kernels representing 
Hankel-Iike integral operators on (xo,oo), one has, using the inequality ||AB||^ < 
I I'^l l«sl 1^1 1«5 foi" Hilbert-Schmidt operators, 

\\Kn-K^\\^^ < ||A*"^(«+3')-.-("+>')%7(«+j,T)|LJ|^i")(v+j)|L 

+ ||.-("+>')U7(«+j,T)|LJ|/("'(v + 3;)-.(^'+>')^A+(v+j^ 



where 



C' = max(||."("+>)U7(« + 3^,T)|U,,sup||^(")(v+j)|L,)C.~-<'. 

n>N 



Then, from ( |2.14[ ), one has that 

2C' 

|Iogdet(/-/:„)-logdet(/-/:ec)| < -n'^^^ 

l-M 



Remembering the representation of the probabilities in the statement (2.13 1 in 



terms of Fredholm determinants establishes Corollary |2.4| □ 

3 An integrable deformation of Gaussian random ensemble with 
external source and 3-component KP 

The connection between the Gaussian random ensemble with external source 
and the multi-component KP hierarchy is explained in lIU and IS. The main ideas 
are sketched in this section. For the multicomponent KP hierarchy, see [30]. 

3.1 Two sets of weights and the p + q-KV hierarchy 

Define two sets of weights 

i/^i(x),...,VA^(x) and 91 ();),..., ^^(j), withx,3;eM, 

and deformed weights depending on time parameters Sa = {sai,Sa2, ■ ■ ■) <0!. < 
q) andtp = {tpi,tp2, ■ ■ ■) (1 < jS < p), denoted by 

(3.1) Ya\x) := \lfa{x)e-^^-''-''' and (p'p{y) := (pp{y)e^^-^'P^y' . 

That is, each weight goes with its own set of times. For each set of positive inte- 



ger m= {nil, • ■ • ) 'fiq), n = {ni,. . rip) with \m\ = \n\, consider the determinant 



=Ia=l™a and |n| = ijj^i "/j. 
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of a moment matrix 7]„„ of size |m| = |«|, composed of pq blocks of sizes m;?iy; the 
moments are taken with regard to a (not necessarily symmetric) inner product (• | •) 



(3.2) TififiiySi, . . . ,Sq,t\, . . . ^tp) 



:= det 



</<m[ 
0<j<»i 



0<j<n| 



0<j<np 



0<j<n„ / 



We now state a non-trivial Theorem involving a relationship between the deter- 
minants of the block moment matrices above, by increasing or decreasing the sizes 
of the blocks by one. Modifying the size rip inn = {ni, ... ,np) by 1 is indicated 
by n n±ep, where ep = (0, ... ,0, 1,0, ... ,0), with 1 at place p. The proof and 
many simple examples can be found in ||6l: 

Theorem 3.1. (Adler, van Moerbeke and Vanhaecke S) Then the block matrices 
imn satisfy the [p + q)-KP hierarchy; to be precise, the functions Tmn satisfy the 
bilinear relation]^ 

(3.3) 

for all m,n,m* ,n* such that \m*\ = + 1 and \m\ = — 1 and all s,t,s* ,t* £ C°° 
and where aa{m) =Y,^,^i{ma' —m*^,) and <yp{n) =J^p,^^^{npi - n*p,). 

Computing the residues in the contour integrals above, the functions with 
\m\ = \n\, satisfy the following PDE's in terms of the Hirota symbol, defined in the 



The integrals are contour integrals along a small circle about oo, with formal Laurent series 
as the integrand. Also, for z e C, we define [z^'] := . . .). For a given polynomial 

p(?i,f2, ■ • ■ )' the Hirota symbol between functions / = /(fi ,f2, •• •) and g = ,?2, ■ • ■) is defined 
:= /'(^'^^'•■•^■^(^ + ^)^(^^-'*')|v=o' ^^^'^ elementary Schur 

polynomials s^, defined by eEi*''^' := Y,k>0^k{f)z'^ for ^ > and S/'(f) = for ^ < 0; moreover, set 
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footnote 11: 



2 



'/3/+lC"/3',l 

,2 



.2 



2 

(3.4) - logTmn = — Sf(5.v„)Tm-e„,,i-t.„ o Tm+e„,n+e„. 



3.2 Gaussian ensemble with external source 

Consider an ensemble of n x n Hermitian matrices with an external source, 
given by a diagonal matrix A = diag(ai, . . . ,a„) and a general potential V{z), with 
density 

P„(M G [M,M + dM]) = — f.-Tr(^W-^^)t/M. 

For a subset £ C M, the following probability can be transformed by the Harish- 
Chandra-Itzykson-Zuber formula, with D := diag(zi , . . . ,z„), 

A„(z):= n fe-^i)' 

l<i<y'<n 

and all distinct a;, 

(3.5) P„(spectrumMc£) = — / e-Tr(y(M)-AM)^^ 

Ai JE" I Ju{n) 

= / A„(z)det[e-^fc)+'''^^] !<,,•<„ n^/z,-, 

with Gj 7^ <3j and the Vandermonde A„(z) = ni<i<;<«(z! The formula re- 

mains valid in the limit, when some a/'s coincide, upon making differences of rows 
and dividing by the appropriate (a, — ay)'s. In the following Proposition, we con- 



sider a general situation, of which ( 3.5 1 with A = diag(a, . . . , a, 0, . . . , 0) is a special 
case, by setting (p+ = e"^ and (jo^ = 1. Consider the Vandermonde determinant 
A„(x,j) := A„(xi, . . . ,Xytj,ji, . . . jjytj)- Then we have the following (see lH): 
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Proposition 3.2. Given an arbitrary potential V{z) and arbitrary functions (p^{z) 
and (p^{z), define (n = k\ +^2) 



(pi,...,p„) -^'-^W (<p+(z),z<p+(z),...,z^'->+(z), 

(p-{z), z(p-{z),...,^^-^(p-{z)) 

We have 



If " 
(3.6) — / A„(z)det(p;(zy))i<,-j<„r[<iz/ 
nl Je" I 



/ Anix,y)At,ix)At,{y)Yl9^ixi)e-''(-'^^UxiYl(p-i^^^^^ 

ki\k2\.JE" 1 I 



det 



0< j<A:i+A:2-l 



0<j <ki+k2-l I 



3.3 Adding extra-variables t, su, and /3 

We add the extra-variables f = (?i,f2, • • •)) ■5' = (^i ,52, • • •)> « = (mi,M2, • • •) and j8 
in the exponentials, as follows {n = +^2), 



(3.7) 



<p-(z) = 



, <p (z) = 



The determinant of the moment matrix ( 3.2 1 with regard to the inner-product (/, g) 



1, q = 2, n[ = k\+k2, mi = ki, m2 = k2, and 



JEf{z)g{z)e ^ l^dz, with p 
(p\{x) = 1, V'^i(3') = e"^^^^ , ^2{y) = 1 is the same as the determinant (3.6 1, with 
the expressions V(z) and (p^{z) as in (3.7 1, and setting ^i, := Si, S2i = Ui, tu = f,-. 
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Therefore by virtue of Theorem 3. 1 the expression below satisfies the 3-KP hier- 
archy, since p + q = 3, namely, 



Zkik2{t,s,u;a,p;E) := det 



E J 0<i<h-l 

0< j <k,+k2-l 



(3.8) 



ki\k2 



Corollary 3.3. The functions x^^f^^it ,s^u) satisfy the identities 
(3.9) — log^ =^^-L- ^ __iog. 



dh Zk,-l,k2 g^logT^,,^/ '^^1- 



l,k2 



log ^^1.^2 
5& log ^^1,^2 



y^^ry. , '^kiM + l _ dhd^i^^^'^k, ,k2 d T:ki,k2 + \ _ dndli2^°^'^ki,k2 



52 



'^^1-^2-1 37§;7logT^,,,t2' ^^^1.^2 



Proof. The bilinear identities ( |3.3| ) imply the PDE's ( |3.4[ ) for 

'^kik2{Us,u) := 'Zkuk2M+k2{^^^^^)i 



(in the notation of ( 3.2 1, setting m = {ki,k2) and n = + ^2) expressed in terms of 
Hirota's symbol, for j = 0, 1,2, . . ., 

(3.11) Sy(5,)T<:i + l,/:, OTi:i_l,A:2 = ~%k2 g^^g^,^^ log^/i:iA:2 

2 



S;(5/)Ti:|,A:2 + l O 



.h-l 



~^lk. 



52 



■logT/t,*:2 
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In particular for j = 0, one finds the following expressions 

, - d^logTk.M _ _%Jc2fi%M-l 

and another set of expressions for 7 = 1. Then taking appropriate ratios of these 
expressions yields the formulae of Corollary [33] □ 

4 Virasoro constraints 

Define the differential operators involving the boundary points of the set E, 

=^m:=£^r+'^, for E = \J[b2i-ub2i]cR, 
'=1 ' (=1 

and the differential operators V_ i and Vq involving differentiation with respect to 
the auxiliary variables ti, Sj, Ui and j8: 

d d ^ ( ■ d ■ d . d 

V-i := -^-215- h> iti^^ \-iSi^ \-iUi- 

dti dsi g V dti^i dSi^i dui-y 

+ki {ti - si) +k2{t\ - ui) + aki 

d d d ^ ( . d . d . d 

Vo := -^-2p^ «^ + L iti^ + iSi^ + lUi — 

dt2 ds2 dsy f^^ \ dti osi dui 

+k\ + kl+kik2. 



Theorem 4.1. The integral Tii^k2{t,s,u;a,l5;E) defined in (3.8) satisfies : 
(4. 1) Tk^ki = VmTi:,A;2 form=- \ and 0, 

with Ym as in (14. 1 P. 



Proof. Let 4j/;j(jc, be the integrand of the integral Zic^k2{t,s,u;a,p;E) defined in 
CT : 

r k[ k2 

z := Zk,k.{t,s,u;a,p;E) = / Ikjk2{x,y)Y\dxiY\dyj- 

■Je" f- J 

To obtain (4J_I, one uses the fundamental Theorem of calculus and the commuta- 
tion relation 

A 5 d 

one uses the fact that for any Vandermonde, 

' d 



Y,^H^u---,xk) = o 
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and 



k 



^ d n(n — \) 
2_^Xi—A{xi,...,Xk) = A{xi,...,Xkj. 



Also one uses the fact that the sum of the Xj and yt derivatives of the integrand 4j 
translates into the tj, Si and M,-derivatives of 4i,/t2> the auxiliary parameters 
t,s,u;a,p were precisely added for this very purpose! So we find: 

r / 3 d \ ''^ 



E" 



^-lihi.k-dWdxiYYdyj 
1 1 



^ r h \ 

Jj,Jkuk2Y\dxiY\dyj\ =V_iT 



and 



h k2 

Vo (4, ,k, ) Yl n 

E" 1 1 



Vo I Ik, M n n "^yj ) = voT, 



establishing Theorem 4.1 Aother way of computing this has appeared in [2]. □ 



We also have the following identities, valid when acting on Tkik2 {t ,s,u;a,^;E)\ 
(4.2) 



d d d d d d d d d d 

(4.3) = , — = , = , — = . 

dsi da^ dti da 5mi ' ds2 5j3 ' dt2 dj5 du2 
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Corollary 4.2. On the locus ^ = {t = s = u = ^ = 0}, the function 
f := logTyt,/t2(f,5,M;o;,j8;£') satisfies the Virasoro constraints : 



(4.4) 



dti 

dt2 
dl 
du2 



da' 



ds2 dp 
du\ 



da 



?o + a^ ) /+ ( k\ + kl + kik2 



(4.5) 



dt\du\ 
dt\du2 
dt2du\ 



da 



da 



+ ^ ]f-2aki 
da 



(4.6) 



dt\ds\ 
dtids2 
dt2dsi 



5-1 



d 

da 
d 

d^ 



f-ki 



%-a 



da 



d 

da 



f 



Proof. Upon dividing by T, equations (4.4 1 are a direct consequence of (4.1 1 and 



(4.3 1, when evaluated on the locus To derive equations (4.5 1 and (4.6 1, we use 
the fact that the boundary operators 



2'- d 

V f,m+l " 



commute with pure time-differential operators. For example, the calculation of the 



first equation in (4.5 1 goes as follows. We know from (|4. 1|) and (4.3 1 that : 



if-- 



dti 



ak\ +L\f + i 



du\ 



+ aki +L2/ + J 
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where the L, are Unear differential operators vanishing on ^ (commuting with 
SBrn) and the are functions vanishing on ^ (/ = 1,2). Therefore: 



-1 + ^ )/ 
da 



+ aki+L2f + i2 



if 



dui \dt\ 
Since 



aki +L\f + \ 



d 
du\ 

d£i 



if 



d^f d 
J(' \dt\dui du\ du\ 



i>2 



i-l 



dsi- 



+ iui 



and l\ = ^1^1 + ^2"i — (^1 + ^2) h, one checks that, along the locus one has 



Li/ = 



and 



yielding the expression for g^J^^^ in (4.5 1. A similar procedure applies to establish 
all the identities above. □ 



5 A PDE for the Gaussian ensemble with external source 

Consider the Gaussian Hermitian random matrix ensemble J^„ with external 
source A, given by the diagonal matrix ( 0.16| ) (set n = ki+ kz) and density 

(5.1) -e-'^^'^M'-AM)^^^ 

Given a disjoint union of intervals E := |J/^=i [^2(-i , b2i] C M, define the algebra of 
differential operators, generated by 



(5.2) ^, = £ftW^ 



2r 

Consider the following probability: 

(5.3) P„(a;£) := P( all eigenvalues G £) = — / e-'^'^^^^^-^^UM, 

where J^i{E) is the set of all Hermitian matrices with all eigenvalues in E. The 
purpose of this section is to prove the following theorem: 
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Theorem 5.1. Tlie log of the probability P„(a;£') satisfies a fourth-order PDF in 
a and in the endpoints bi, ...,b2r of the set E, with quartic non-linearity: 



(5.4) 



F^G -^F G+ 



= 0, 



or what is the same 



(5.5) 



where 



det 



da 



( 


G+ 


m-iF+ 


-F+ 





\ 




-G 


m^iF- 


-F 







c 


?_iG+ 







-F^ 




\ - 




SB\F- 





-F 




logi 




F- = 


-^-i ( 




+ - 
< 



0, 



<9a 



logP„-A:2 



Hi 



d 

4— logP„+4aA:i + 
da 



4kik2 



a 

d \ d 
a— + 1 1 — logi 



da 



?o-a 



da 



da 

-la^-x 



a 



2G 



%-a- 



-1 



da 
d 



log I 



da ) \ da 

{lif.F^],,_^^{Ht.F^],„,. 



log I 



Remark 5.2. It is not surprising that the PDE ( |5.4[ ) has exactly f/ie same form as 
the PDE derived in [4] and |5|, associated to the Gaussian Unitary Ensemble with 
an external source, in the case where the source matrix admits two eigenvalues 
of opposite signs. The only difference is that the expressions for the functions 
F^ ^H^ and G^ obtained here, differ from those in The reason is that cor- 
responding T-functions satisfy the same integrable equation (the 3-KP hierarchy. 



as in Section 3.3), whereas the Virasoro constraints leading to (5.4i are different. 



In particular, there is no more involution relating the variables Sk and Uk- 



Proof of Theorem 5. 1 Remember: 
P„(a,£') := P„(specM C E) 



Tkik2{t,s,u;a,[5;E) 
Tk,k2{t,s,u;a,p;R) 



^={t=s=u=p=0} 



The denominator, that is the integral ( 3.8 1 over the whole range, can be expressed in 
terms of moments, which contain standard Gaussian integrals (the reader is referred 
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to the Appendix 1 in lU), leading to an exact evaluation, with Cyti/tj a constant, 
depending on /c 1,^:2 only: 



(5.6) 



Zkuk2it,s,u;a,l5;l 



Consequently (with Qjj;, a constant, depending on ^1,^2 only): 



(5.7) \og¥n{a,E)=\ogT:kuh{t,s,u;a,^;E) ^ - —a^ -kik2\og{a) -Ck^ko- 



Then, we turn our attention to the numerator T/:,i2(?,'5',M;oi,j8;£'), more briefly 



noted as 'Ck^ki- On the locus by (3.9l, and (4.4i, one finds: 
(5.8) ^^"^"^-'^^ = ^log f = -^-ilog ( 1 +2a 



logT^jl, dti \T:ki-\,k2j 



dtidsi ""to '"1 
dttds 



^ki-l,k2 



log 



dtidsi 



j-logTk,^k2 da \Zk 



log 



-\,k2 



In (5.8 1, acting with d/da on the first expression, then acting with on the 
second expression, and adding the two, yields: 



This identity can conveniently be expressed as Wronskians: 



0. 



(5.9) <!=^^_i^logT^,^, 



T 



3_ 
da 



where 
(5.10) 



7+ .- 



■ log % 



dt\dsi 
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In this way, by virtue of (4.6 1 and (5.7 1 one obtains, along the locus Jf, explicit 
expressions for F+, and H2 that are free of partials in j3; namely: 



(5.11) F+ 



Hi 



da 



log Til -^1 = 



da 



4(^)logT,,,=4(^)logF„+4a..+4^ 



%-a- l-2a^-i 

da 



%-a^ \-2aSS-i 

da 



d 

da 



da 



logTiii2+4a^i 



logP„. 



Subsequently, one repeats exactly the same operations for the Virasoro and KP- 
identities involving the and Uk variables. From (3.10l and (4.4i, one finds on 



(5.12) 



dt 



;a^logTi,l2 d 



^logT,,.,, dt 



d 

dt{au{ 



-logf^^U-^^ilog 



du 



l0£ 



'^ki,k2+i 
'^h,k2-\ 

da 



So, in (5. 12 1, acting on the first equation with 
equation with 1 , then subtracting the two yields: 



^-1 + ^ 



'^<:i,/t2 + l 



'^h,k2-\ . 

and on the second 



dtidiii 



dt 



52 



+ 



d_ ( log ^^1.^2 ' 



0, 



or equivalently (remember the brackets are Wronskians) : 
(5.13) 



H 



1 ' 



+ {H. 



'2 ' " 



3 



in terms of the functions 

d^ 



(5.14) 



F 



dtidui 
d^ 



^2 \ g 

-^TT^ 373~ ) log ^kt M - 2-^- 1 log Tkik2 



dtjdui 



^OgZk,^k2- 
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Using the Virasoro constraints ( |4.5| ), as well as ( |5.7| ), we obtain explicit formulae 
for F , and H2 (which do not contain partials in j8): 

d . 



^-1 + 



da 
d 

da 



HZ 



(5.15) 



d 

da 
d 

da 
d 

da 
d 

da 



da 



%-a— + \ 
da 



%-a 

aa 
d 

%-a- 1 



logP„-^2 



fVogP„-4^ 
da a 



+ — )logP„. 
da ' 



-4aki 



Equations ( |5.9[ ) and ( |5.13| ) form a linear system in 

and 



dlogZkik2 



dp 



2 dlogTk.ko 



X ' dp 

which can be solved to extract the quantities: 

dlogrk,,k2 ^ G F+ + G+F 

dp y -F-{^-iF+) + F+{^-iF-) 
,2 d\ogXk,,k,\ _ G-(^_if+)+G+(^_iF-) 



(5.16) 



(5.17) 



dp 



-F {^^iF+) +F+ i^-iF-) ■ 



Finally, subtracting the second relation from 1 of the first equation, establishes 

To pro 

^OgTkik, 



the expected PDE ( 5.4 1 and Theorem 5.1 To prove the second equation ( 5.5 1, set 

. d 



X :-- 



dp 



then the matrix in (5.5 1 annihilates the column (l,X,^_iX,=^ijX) , and thus the 



determinant vanishes, concluding the proof of Theorem 5.1 



□ 



6 A PDE for the transition probability of the r-Airy process 

Proof of Theorem 0.4 Remember from section [T] the identity 



Um 



j{0,p7«72) 
Br 



all Xi 



1 



l+e-2V«' 



2cosh(T/?iV3) / 
limP„ (p^e^l""';l^ + 



„l/6 
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which for < p < 1 leads by Theorem |0.1| to a phase transition at p = 1 , for which 
the expression above reduces to P(sup.^^(T) < x), according to (0.7 ). The above 
scaling suggests the choice z = n-^^^ as small parameter and considering the map 

(t,x) I — > {cc,b), given by 



a = py/ne' = ^e^' 



b = 2^/n + 



,1/6 



■ +XZ 



with inverse map {a,b) i — > given by 



1 , / az" 
T = ^ lo£ 



Z z 



4 ■ 



Setting 
(6.1) 



Q{t,x) := logP„ ( p^/ne'/"''';2Vn + 



„l/6 



logP„( ^^'^^;^ + (-oo,x)^), 



yields, via the inverse map, 
(6.2) 



1 , ,az\ b 2 
Q{ ^log 



Z2 ^> ^'z z4 



logP„(a,(-oo,Z7)). 



From Corollary |2.4| (section[2]) it follows that for z — > 0, 
(6.3) Q{z,x) = Q{T,x) + 0{z^), 

with independent of z. 



orem 



Si nce we have shown that logP„(a, (— oo^^)) satisfies the PDE (5.4 1 of The- 
with = ^ and = > it follows that, to get a PDE for the 



5.1 



limiting case, we just need to estimate: 



(6.4) 



F+G- +F-G 



b ^xz+^ 



To do this, the various expre ssio ns in the bracket must be computed in terms of 
the function Q{'C,x). By (6.1 1 and (6.2), one immediately gets: 



3logP„ ^ 1 fdQ 
az^ V dx 



—e 



dQ\ 5logP„ IfdQ 



da az^\dz J p" \dz J' db z\dx )' 

Hereafter, to shorten notation, we will write the partials as 



dxdx 



Q, Gt« 



dT^dx 



Q, etc. 
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Without taking a limit yet, but expanding asymptotically the expressions in powers 
of z, we find 



F+g§^F~-F' 



Qrxx 




pz' 


pz^ 


Qtxxx 




pz^ 


pz 



pz^ 



+ - 



pz" 



In o rder t o compute the expansions of (which are respectively defined in (5.9 
and (|5.13| ) and of =^_iG^, we need the asymptotics of Hf, as defined in (5.11 
and (pTra. In the end, one finds 



2Qxxx 


2Qtx + r{p- 


2 


' Qxx[{p- 




\ -2rp{p- 


2Q^xxx 
p2 


2Qrx + ril- 


2 


' Qrx(^{p- 


~~p^ 





1 

AO 



'7 1 



+4+0(4: 

v z? 



\-2rp{p-l ){Qrrxx + TQr,^,,) J 



where ^ and are given by the following expressions and where and 
denote further expressions in the derivatives of Q, 

Sr = 2?Qy^ + 2r(e„^(2,,, + t2 - x) + 2Te„., + 2Q^, + Q„^) - 2Q\^{Q,„ + 2) 

-2e«(TGT« + QzTX - Qtxx{.Qxx+^1^ +x))+2QrxxQT + QzzQztx + (p - 1)^1 , 



=2.r Qxxxx ~\~ '^f{Q'Cxx.x{Qxx ~^ ^ —x)-\- 2t2ttxv + Qtxx + 2tttx) 

~^QtxQxxxx ~ '^Qtxi^Q'crxx + Gtttjc + ^Qtxx ~ QrxxxiQxx + 4t + x) ) 

Consequently, using the two leading orders for the expressions above, one obtains 
for small z- 




F+SS^iQ- +F-ii8^iG 

- (f+g-+f-g 



F+^l^F- -F-^^^F+ 



4(P-1) 



2rp- 



dxdx) \ dx^dx' dxdx^ 



m)+(p-i)^{Q) 



O 
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where S'{Q) is given by 



^(G) = 2 r 



(6.6) 




+2 r 



+ 



drdx^' dr \ dr \ dr dx 



+ 



dxdx^ 



dl V dx 



where ^{Q) is a similar expression, that will not be needed, and where the bracket 
is a Wronskian associated to the "space" operator d /dx. 



Hence, /or < p < 1, taking the limit z — > above and using (6.3 1, yields the 
equation 



2rp 



drd} 



dx^dx ' dxdx^ 







which is trivially satisfied; indeed, from other considerations we know that is 
the Tracy- Widom distribution, which of course is T-independent. 

However, in the critical case p = 1, the leading term has order 1 /z'^, with co- 
efficient <a{Q{T,x;z)), as in (6.6 1. Then taking a limit in (6.5 1 when z — > and 
using ( |6.3[ ), forces upon us the equation 

lim^((2(T,x;z)) = ^(e(T,x))=0, 

which an easy computation shows can be written as 
(6.7) 



dzdx^' 



.2 53 2 



+ r 



'dQ 



+ 



1 



1 / d^Q 

2 [drd^ 



dr^^dx^^^ dx ) dTdx^ ^ \ dtdx ' j- ^ 
d^Q d^Q d^Q I ( d^aV d^Q I f d-Q djtQ) 1 

dl^ dldx dx^ 



and further rewritten as equation (0.12 1, ending the proof of Theorem 0.4. 



□ 



Remark 6.1. This "phase transition" at p = 1 is completely analogous to the results 
found in HI and ||29J for small rank perturbations of random Hermitian matrices. 

7 Remote past asymptotics 

The aim of this section is to study the behavior of the r-Airy process s^^''\t) 
for t —>■ — oo, as stated in Theorem 0.7 In this section T will be systematically 
replaced by t. This theorem will be rephrased as Theorem |7.1[ which includes 
some additional details. 
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Theorem 7.1. The log of the probability of the r-Airy process Q{t,x) = 
log¥(sup £/^''\t) < x) admits an asymptotic expansion, as t ^ — oo, having the 
following form 

Q{t,x) = Qo{x) + jQi{x) + -^Q2{x) + ..., 
for the initial condition 

poo 

lim Q{t,x) := Qo{x) := logP(sup^/(f) <x) = - (a -x)g^{a)da, 

t^-^ Jx 

and where 

(7.1) P r ■ 1? 

Qi = rQ'o, Q2 = -Go, 23 = ^e;," + 3^20, 24 = ^ejj' + y^Go + 3^20- 

5 3 T 3 

25 = ^Go + j^^Qo + 1^20 + 5 

26 = ^2;i' + ^xeii> J^e;$' + ^[,^^ + A(.,2g;;+l3(,+,^)e.■ 



F"-! 1 



nl^^ 3 12(;j-4)! 



1=2 



/or wme constant C(, and with 

noa poo poo 

^5 ■■= x^Q'o + 4^Go + Go + 10 / Qq-6 dy duQ'^. 

Jx Jx Jy 

Also 

P(sup=c/W(0 < x) = P (^sup^(0 < (x+ J)(l + ^) + ^) 
r/ie mean and variance of the right edge of the process behave as 

E(sup^W(0) = E(sup^W(o) (i - ^) - 7 - 

var(sup^/W(0) = var(supi;/(°)(0) (l"^) +'^(^)- 

Remember the r-Airy kernel 

k'^''\u,v) = I dwA':{u + w\t)Al:{v + w\t), 
Jo 



as in (0.9 1, where Af{u;x) is given by (0.8 1, where C is a contour running from 
oogS'TT/b j-Q ^gi7c/6^ gy^j^ jj^^j jjes above the contour. In this section one lets t — > 
— oo, which, of course, implies that —it will remain above the contour C and thus 
this limit is compatible with the contour just mentioned. Letting t +oo would 



40 M. ADLER, J. DELEPINE, AND P. VAN MOERBEKE 

require a drastic change of tlie functions Af. In this section the subscript t will 
occasionally be omitted from the r-Airy kernel Kf'\u,v). Note that iC'/"'(M,v) = 
A'(")(m,v) is the Airy kernel, which is independent of t, 



A{u)A'{v)-A'{u)A{v) 



u — v 



A(w + u)A{w + v)dw, 



where 



A{u) =Aq{u) = e3 



In 



is the Airy function, satisfying the ordinary differential equation A"{u) = uA{u) 
and behaving asymptotically as 



(7.2) 



A(x) 



_2 3/2 

e i 



at 



(;c3/2) 



T + . . .), as ■ 



A^{u)du = {A'^-AA"){x) 
j^A^{u)du = ((a2)'+x(a'2-AA")) (x). 



The ODE and \mix^ooA{x) = imply the following formulae, upon differentiation 
by X, 



(7.3) 



Also remember the Tracy-Widom distribution 11261 

det(^/ =exp(^-J^ {a-x)g^{a)da^ , 

where g{cc) is the Hastings-MacLeod solution (0. 14 1 of Painleve II. The following 
shorthand notation will be used for integers £> I, 



O(A^) = O 



( 



\ 



e 3 



7 3 
2^2 



20rxV4 



for X 



whatever be the power ^ G M. 
Lemma 7.2. Given 

eoW = logdet (7-^(0))^ = 
one checks 

poo 

Q'o{x) = / g^{u)du=A^{x)-A"{x)A{x) + 0{A 



a-x)g^{a)da, 



Q^^\x) = -{g'f"-'){x) = -{A'r-'^{x)+0{A')Jorn>2. 



,2^(n-2) 
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Proof. The estimates follow from ( |7.3| ) and an improved version of the estimate 
( |0.14| ) by Hastings-McLeod L18J, namely: 

/ _ 4^3/2 \ 

g{x) = K{x) + O I J for ^ / oo, 

which is to be interpreted as a genuine asymptotic formula; i.e., it can be both 
integrated and differentiated. □ 

For future use, one needs the following estimates for the Airy function: 
Lemma 7.3. For x ^ oo, one has the estimates 

[A\u)du = (a'^-A"A)W = =0{A^). 

Proof. Upon using integration by parts and upon substituting the asymptotic for- 
mula ( |7.2| ) for the Airy function, one computes for instance, 

100 ~/ 



1 

e"3" . ^ c': , 1 re^i" , ^ c' 



Further terms in the expansion can be obtained by differentiation by parts and sim- 



ilarly for the first expression, thus ending the proof of Lemma 7.3 □ 



Given a kernel F{y,u) acting on L^{E) with £" C M and a bounded continuous 
real function / on define the norms 



F||i = sup / \F{y,u)\du and ||/||oc = sup 1/(3;) |. 

yeE JE yeE 

If / is a function of several variables, the sup is taken over all variables. Then 

FG\\\ < sup dz. du\F{y,u)G{u,z)\ 

yeE JE JE 

< sup / du\F{y,u)\&up I dz\G{u,z)\ = \\F\\\\\G\ 

yeEJE ueE-JE 



yeE 

HenceEi] 

00 00 I I I I 

l^lli 



1 1 A 



1 1 



^2 Obviously the estimate below requires ||-F||i < 1. In the application, this is achieved by re- 
stricting the domain of the operator to the interval (x, <x) for sufficiently large x. 
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and 

||F/||oo = sup / \F{y,u)f{u)\du < sup|/(M)|sup / \F{y,u)\du = \ \F\ 

yeE JE ueE yeE JE 

Given the kernel A'^"' , define tlie resolvent kernel R by 

(7.4) I+R:={I-K^°^yK 
Then readily 

(7.5) R-K^^^ = i^(«)2(; + ^(0) +^(0)2 ^ ^ _ ^ (/-i^«'))-li^{0)2. 



Lemma 7.4. One has the following estimates 

= 0(a2), ||ir(*')||oc = 0(A2), and \\R- K^^^\\^ = 0{A^). 



Proof. Using the fact that A{u) > is monotonically decreasing for m > 0, and 
setting E = {x, o°) , 



= sup / dv / dwA{u + w)A{v + w) 
ue{x,oo)Jx Jo 

poo pOQ 

< sup / dvA{v) / dwA{u + w) 

ue(x,oo)Jx Jo 

< ^ dvA{v)^ =0{A^) 



by Lemma 7.3 while 



sup |A'(m,v)| = sup 

M,l'£(,V,oo) H,ve(.V,oo) 



A{u + w)A{v + w)dw 



< / A{x + w)^dw 
Jo 

poo 

= / A{w)^dw = 0{A^ 

Jx 



also by Lemma 



7.3 
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Considering the function R{-,v) for fixed v, one has, using 



|«{-,v)-ir<»)(.,v)|U = lUi-Ki'Y'iK'o^H-MW- 
<||(;-x<°')-'||,||;fi°i^(.,v.) 



Hence, with ||A'||oo := sup |/r(M,v)|, 

u.ve£ 



|/?-^W|U:= sup |/?(.,v)-/^(°)(.,v)|<-^^^||^(«)||. = 0(A'' 



ending the proof of Lemma 7.4 □ 



The next point is to get an asymptotic expansion for the Fredholm determinant 
det{I — Kf '^^'')x,oo in t, using the asymptotic expansion of the kernel Kj'\u,v) in t 
(Lemma |2!2l ); this leads to the next Lemma: 

Lemma 7.5. The following probability for the r-Airy process has an asymptotic 
expansion in l/t of the following form 



Q{t,x)=logF{sup^^'-\t)<x) = logdet(l-K^'-A =2oW + £^, 

^ n=l ^ 



where 



An 



(7.6) Q„{x) = £ r"-^'Qn^„^2i{x) = - TvK^ + 0{A^) , 



i=0 
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where Qq = logP(supi2/(f) < x) and where the Qi — > and have all their deriva- 
tives Oforx oo. Moreover, for X o°, one has 

Qi{x) = r{A^-AA'' + 0{A^)) 
Q2{x) = -^{A^ + 0{A')) 

Qsix) = --{{A^y + 0{A'))--{{A^y + 3 A'\u)du + 0{A'))+0{A') 



!2«W = -^((a2)(«-2)+0(a4)) 
nl 

•1-2 /T,ri-] w-1 \ ("-4) 



(n-3)! V 24 ' ' 2 
+ £ r"-^%,„^2i{A) + 0{A'), 



1=2 



with r„ „_2((A) = quadratic polynomial of A and A', with coefficients depending 
onx -\- (quadratic polynomial of A and A', with coefficients depending on x). 



Proof. We shall always operate in L^(x,oo), so that occasionally the x will be sup- 
pressed. Then, using the asymptotics for the kernel Kf""^ as in Lemma 2.2 one has 
the following: 



t t ;^ 1 ^ 



I>1 



with (the resolvent operator 7? of the Airy kernel is defined in (7.4 1) 
(7.7) Li = = {I + R)Kj''K 

Using log(l — z) = — z — y — y — . . ., one finds 

= logdet [I-K^''^)^ = Tr\og{I-Kj'-'>) = Trlog(/ + £ 

where 



,=1 ^' 



(7.8) Gi = -TrLi, (22 = -Tr(L2 + ^), & = - Tr^La +L1L2 + ^ 

24 = -Tr(L4 + ^lI + L]L2 + UL3 + ^Lt) , . . . . 

More generally, the 2„'s are weight-homogeneous polynomials of degree n in the 
Li, with weight(L,) = /, having the form below, which can further be expressed in 
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terms of the K- and R, using expression (7.7 1 for the L;, 
(V.9) Q„ = -TrL„ + TrP„(Li,L2,...,L„_i) 

= - TrKj[^ - Tr{RKj[^ ) + Tr /'„ (Lj , L2 , . . . , L„_ i ) 

= - Trd''^ + TrSniK^'^ , . . . ^K^'l^R), 

with Pyt and Sk polynomials of non-commutative variables with no linear or inde- 
pendent terms, but with quadratic terms and higher. From Lemma 2.2 the kernels 
K-'^ and hence the L,- are polynomials in r of degree / having no constant terms; 
hence the 2„'s, by their weight-homogeneity, are polynomials of degree n in r, 
having no constant term and so (here one must indicate the r-dependence) 

r' 



(7.10) 



!2.(?,x) := e(?,x) = £ (- j £ 



The claim is that only the terms Qn,n-2i{x) appear. Observe from (0.9 1 and (0.8 1, 
that k'^/'^ (m, v) = ^/'^ (v, u) , and thus 



Qr{t,x) = d&t{I-K^''\u,v)), = det{I-Kr{v,u)), = det{I - K'_/> {v,u)) 

= Q-r{-t,x). 



So from (7.10), one has 



1 



Qr{t,x) = -{Qr{t,x)+Q-r{-t,x)), 

implying that only the even terms appear in the sum in ( |7.10 l, thus proving the 
statement ( |7.6| ). 

Next we now proceed to estimate the two traces above: 
• At first 



TrS,{Ki'-\...,Kl''^„R) = 0{A'), 

which we now illustrate on a typical example, like Tr{RK\ ) . Note A (m) > for 
M > 0, and so by Lemma [74] 



('•) 



I TtRK 



1 I 



Tr( {R - K^^^ )k[''^ ) I + I TriK^^^Kl'^ ) | 



< (WR-K 



\l.lx 

< (0(a4) + 0(a2)) jj^ A(u)A{v)dudv 



< 0{A^ 

< 0{Ay, 



A{u)du 



the last estimate follows from Lemma 7.3 More generally, the trace of a monomial 
of degree i has order 0{A^^). 
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(r) (r) 

i Then we evaluate Tr^,, ; in order to do so, it suffices to evaluate the kernels K„ 



of Lemma 2.2 along the diagonal, and to notice that a skew-symmetric operator 
vanishes on the diagonal. Since the domain of the operator is unbounded, one needs 

(r) 

to consider K), Z(v,„)' trace is then obtained by integrating on the diagonal and 
by taking the limit m — > oo, upon using the decay of the kernel at oo. Therefore, on 
the diagonal, Kn \u,v) is a polynomial of degree n in r, skipping every other term. 



kI[\u,u) 



r"-^ f3n-\f n-\f , a ("-3) 



(?i-3)! V 24 



+ . 



That the Qi and that Qo = logP(sup < x) has all their derivatives 
for X ^ oo follows from the statement on Tn.n-niJ^ (which are the coefficients of 

) and the asymptotics of 

□ 



r" appearing in Qn{x) as in the statement of Lemma 



7.5 



the Airy function, from which Lemma 7.5 follows. 



Proof of Theorem 7.1 From Sectionj6j we know that 2 (f,x) = log P(sup (f) < 
x) satisfies the non-linear PDE (0.12 1; it is more convenient here to use version 
\6.1\ of the equation. Also remember Q{t,x) 0, when x ^ oo. Then, assigning 
weight= 1 to both variables t and r, one readily checks that the PDE (6.7 1 can be 
graded as follows: 



d'Q 2 




dtdx^''' 


5x3 






\ dtdx^ ' 





(weight= 1) 
dQ, d\xQ) 



dt^^^ dx' 



dtdj 



+ 2 



\dtdx' 



-other terms of weight < —1. 



dx^ 



(weight= —1) 



Since, by Lemma 7.5 the solution has the following general form 



Q t Q \t / \ r r 



it follows that one can compute inductively all the Qnn {x) and then inductively all 
the Qn,n-2{x) and so the Qn will be as announced in Theorem |7.1[ 

Setting this solution in the PDE above, yields a series of descending weights, 
namely = Wi + W_i + W-3 + which holds for t -00, all x G M and all 
integers r > 0; this implies Wi = W-i = W-s = ... = 0; one then checks the explicit 
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expressions 



^ = -i^i:{("-^)e;:-,„-;,e-} 

n=l ' 7=0 



^1 n— 1 
^ ^n+l ^ 



n=l 



{n-j)Q:_j^„_j, 



-2 £ i{Q'u,Qkk} 

e+k=j-\ 



I 



Since this holds for all t 



\ — oo and r > integer, one must have for all x 



J 



(7.11) 



I(«-7){e:-,,„-y,e^} = o, 



1,2,... 



with goo = Go(-x)- 

• For n = \, this is {Q'/j , Qq'} = 0, leading to 2i i = cqQ'^ + ax + ^. Considering 
the asymptotics for x oo and using Lemmas 7.2 and 7.5 a and j3 must = 0, 
leading to the equation 

= 211- coGo = (a'^ -AA")[\ - Co) + 0{A^) = (1 - co)0(a2) + ^(A^), 
implying co = 1 and so = g^. 

• For n = 2, the equation, with the previous data introduced, reads 

= { 2G^'2 , } + { eT 1 , G'/'i } = { , - 2e^'2 + eiJ' } , 

which upon solving leads to Q22 = ^Qo+ciQ'q + a'x + j8'. For the same reason as 
before a' = = 0. Then again using Lemma 7.2 and Lemma 7.5 one finds for 

X ^ 00, 



= Q22-\Qo-ciQ'o 



-f + f -ci(a'2-A"A) + 0(a3) 
ciO{A^) + 0{A^), 



implying ci = and thus Q22 = \Q'^- 

• By induction, assume Qu = j\ Qq\ for < / < « — 1. Then substituting this iden- 
tity into equation (7.11 1 and setting Q„n = ;^2o"' ~^^n leads to pairwise cancella- 
tions in equation (7.11 1 with only one remaining contribution {2o",/?"} = 0, with 
solution Rn = c„2q + a"x + j3", where a" = j3" = 0, and thus, by the asymptotics 
of Lemmas |7.2| and |7.5| 



= (2« 



1 



1 



.(A2)(«-2) 
CnO{A^)+0{A') 



1 



2\ , rttAS 



2^{n-2) 



-Cn{A'^ -AA") + 0{A^ 
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leading to c„ = 0, completing the proof that 2„„ = ^2o"^ for all n = 1,2, .... This 
proves the form of the leading term (coefficient of r") in formulae ( |7.1| ) for the 
2„'s. Since from Lemma 2/2 from the form ( |7.9[ ) of the Q„ and the fact that the 
coefficients Kjf\u,v) in the expansion are divisible by r, the Qi themselves 

are divisible by r. Since they skip every other degree in r, this shows the formulae 
for Qi and Q2; in particular Q20 = 0. 

Setting this information 2„„ = ^q[)^ into the equation W_i = and noticing 
that the following term vanishes automatically, Y,£+k=j-i ^{Q'u^ Q'tk} ~ ^' ^"^^ f^^ds 
for n>3, 



(7.12) 



n-l ( An-j+'i) 



7=0 



n-3 

I 

7=0 



(7+3) 



/ 2L;-l,«-y-3 - (« - j)Qn-j.n-j-2 



\ +(«-7"-3)! 



("-7-2) 



+ 2(«-7-3)e, 



(^«-7-3) 



i For n = 3, by using the fact that Q20 = 0, the equation reads {Qq, {xQ'q — 323i)"} 
= yielding 23i = (5 +C3) + a"'x + j3'" with a'" = {5'" = 0. Thus, using 



Lemmas 7.5 and|7.2[ and the identity (|7.3[), 



= 231 -(3 +4] Go 
= -i ((a2)' + 3j^V(«)^«) - +c'3) (a'2-A4") + 0(A^) 

= - (^£A^{u)du+^(^{A^y + x{A'--AA")')^ -c'^{A^-AA") + 0{A^) 
= -C3(a'2 - AA") + 0{A^) = C30(a2) + O(A^), using dTJ] ) 



yielding C3 = 0, and thus Q31 = |2o- 



• For n = 4, using the formula for 231, the equation reads 



= jeo", (-4242+^ (xeo)'+xe[;+2e'o) | 



AIRY PROCESS WITH OUTLIERS 49 

with solution Q^2 = + (-^+^4)20 + a"x + j8"' and thus a"' = ^" = 0, and 
by the same Lemmas 7.5 and [7^ 

= Q42-\xQl-{j^+c'}j Q'o 
= jl(AVAA") + 3AA"- iA^-AA") + 0{A') 



= -c'^{A^-AA") + 0{A')=c'^0{A^) + 0{A'), 

implying C4 = 0. 

• Using induction, assume 



3(/-3)! 

holds for / = 3, — 1. Then setting 



12(/-4)!^" 



('•-3) 



(7.13) 



Qn.n—2 



3(?i-3)! 



7 



n{n-4) 



into equation (7.12 1 gives the simple equation for /?„, namely {2o",^"} = 0, and so 
Rn = c',2o- Then, rewriting (7.13 1 and using the asymptotics for Qn,n-2 (Lemma 



7.51), and for the derivatives ^ (Lemma 7.2|l, and using the ODE for the Airy 



function xA = A" , we get 

= Qn,n-2 



1 



12(?i-3) 



6{n - 1)a'2 - {3n - 1)(a'2 +AA") 
12(«-3)! V +4AA" -(3?i- 5) (a'2-AA") 

-c',XA^-AA") + 0{A^) 

_c;(a'2 -AA") + 0{A') = c'MA^) + 0{A'), 



(,1-4) 



implying c'„ = 0. Thus the Q„'s are as announced in Theorem |7.l[ namely 



^" + 3 («-3)! ' 12 («-4)! 



lower degree 
terms in r 



In the same fashion we compute and ^52; for example, setting 



into the equation ( |6.7| ), one finds the following differential equation for Q51, namely 
{ Q'^,5Q';, - SxQ'^ + 4Q'^ } + 2Q'^\Q'^ + x)=0, 
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1 



which upon solving leads to 

_ 1 / {x^ + C5)Q'o + 4xQo + Q'^ 

with a constant cj, which has been shown by Aminul Huq (private communication, 
2008) to be 0. Similarly one finds a differential equation for Q(,2 and upon solving 
one finds, for some integration constant ce, 

262 = I (^5 + ^(^'Go + 13(^ + C6)G[,)) . 

Assembling all the pieces, one notices that two Taylor series in Qq and Q'q make 
their appearance in the 1/f-expansion of Q{t,x), leading to shifts in the argument 
of Qo{x) up to order 5: 



(7.14) 



Q{t,x) 



Qi{x) 



+ 



f ^ /In 



xr 7r 
3?3 +T2?4 



7r^ 



xr 



Exponentiating ( 7. 14 1, remembering that e2o(>') = P (sup £/{t) <y), taking a deriv- 

dy 



ative j-P(supi</(f) < y) and setting Pq{x) := P(sup£/(f) < x) yields 



dx 



P(sup^('''(0 <x) 



1 + 



'ifiJdy 



'(sup =6/(0 <y) 



1 



y=(x+i)(l + - 



(7.15) + + + 0(^). 

For the moments with regard to the density ^P(sup^'''''(f) < x), 



°° d 

x^— P(supi^W(0 <x)dx 

— oo 



one reads off from (7.15 1 the following expansion for pil\t) in terms of the t- 

x^PQ{x)dx, 



independent moments 



x''—F{sup£/(t) < x)dx ■■ 

— oo CIX 
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namely, 



H^/\t) = Air 



+ 



ir^iie-l){i-2){e-3)^e-4 + 



'-\){e-2)^f,-3-2e^i) 

14)At£-i) + 0(^). 



In particular, the mean and second moment behave as 

„2 



and 



IrlXi r 
M2-^ + - 



/I 



3t^ 



+ 



6?4 



+ 



Hence the variance of the right edge of the process behaves as 



var(sup=«/M(0) = (Mi^^-Atr^')(0 = ( 1 



2r\ (\ 



var 



(supi;/(°)(0) ( 1 



2r\ f\ 



ending the proof of Theorem 7.1 



□ 



8 The r-Airy process, an interpolation between the Airy and Pearcey 

processes 

Consider n non-intersecting Brownian motions on M, with < < 1 and b <a: 



all Xj{t) G £ for 1 < 7 < 



all xj{0) = 
pn paths end up at a at f = 1 
(1 — p)n paths end up at at ? = 1 



It is intuitive that, when n —>■ oo^ the mean density of Brownian particles has its 
support on one interval for f ~ and on two intervals for f ~ 1 , so that a bifurcation 
appears for some intermediate time to, where one interval splits into two intervals. 
At this point the boundary of the support of the mean density in (x,?)-space has 
a cusp; see Figure 0.3. The Pearcey process describes this cloud of particles near 
the point of bifurcation, with time and space stretched in such a way that the outer 
particles appear infinitely far and such that the time horizon ? = 1 is at infinity. In 
[y it is shown that the same Pearcey process appears in the neighborhood of this 
cusp, independently of the target points a and b, and the number np of paths forced 
to a, showing "universality" of the Pearcey process. It is convenient to introduce 
the parametrization of p, 

(8.1) p = — ^-^ with < <7 < oo. 

Setting for simplicity b = 0, one has the following: 



52 M. ADLER, J. DELEPINE, AND P. VAN MOERBEKE 

Proposition 8.1. 1 1 1 For n ^ oo, the cloud of Brownian particles lies within a 
region, having a cusp at location {xQ^/nfy), with 

{2q-r)a , ^ q2_q+\ y' 

(8.2) xq = — — fo, tQ=\\ + 2a . 

q+l V / 

Moreover, the following Brownian motion probability tends to the probability for 
the Pearcey process ^(f); 

lim P^°/^"^ (^11 Xj (^0 + (coM )^ G XQn^/^ + cqAt + co/l Jy^^ 
(8.3) =P^(^(T)n£ = 0), 



with constants expressed in terms of (8.1) and (8.2 1, 

fq^-q + l\^'^ ^ I toil -to) ^ 



A:=gi/2(i_^)_^-i/2^. 



a a 



The r-Airy process is an interpolation between the Pearcey process and the Airy 
process, which can easily be described by looking at Figure 0.3: 

Theorem 8.2. When — > and for n very large, such that pn equals a fixed integer 
r>0, the tip {xo\/n,to) of the cusp (as given by (8.2)) moves towards the right hand 
boundary of the picture, and, in particular, to the tangency point of the line through 
{a^/n,\) tangent to the curve y = y^2nt{\ — t): 

{xoVn,to) ( I ^2^2 T+W) ^ ^'^^^^y = \/2?i?(l -f)}- 

Also the Pearcey process near the cusp tends to the r-Airy process in the neighbor- 
hood of the point of tangency above. 



Proof. Indeed, letting p ^ 0, or what is the same from ( |8.1[ ), letting q ^ one 
that the cusp is located 

xo\/n = 2ato\/n and to 



sees from formula ( 8.2 1 that the cusp is located at the point 

1 



l+2a2- 

This implies that the point {y,t) = (xo^/njo) belongs to the curve y = y^2nt{l — t) 
and that 

^ ll-to /- 



establishing the first part of Theorem 8.2 That the Pearcey process tends to the 
r-Airy process will be done elsewhere. □ 
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9 Appendix 



The purpose of this appendix is to show that the first few Qi of Theorem 7.1 



can 



be obtained, with tears, by functional analytical methods, in the style of Widom 
|[3n . The proof requires many intricate identities involving the kernels K^'''^ , some 
of which can be found in Tracy -Widom Il26l . This section should convince the 
reader of the usefulness of the PDE's in computing the asymptotics for t — oo. 

Remember the L, = {I + Rx)Kj'''^ from (7.7 1, where we now indicate the explicit 
dependency of the resolvent = K^^\l — A''-"')^^ on x, si nce a ll operators act on 
L^{x,oo). Then Q{t,x) has an expansion, with 2,'s given in (7.8 1, 



Q = mog{i-Kj'-^) = I f = Triog(/ + £ 

(=0 ' i=l 

Throughout this section, we shall be using the inner-product 



Proposition 9.1. 



{/,§)■= / ^{.x,oo){u)f{u)g{u)du. 



2 3 



Lemma 9.2. 



'IrL'l = {-rr{{I + RM,A)". 

Proof. Indeed, 
Tr((/ + /?,)^i)" 

= (-r)" / du:...dun{{{I + Rx)A{ui))A{u2)){{{I + Rx)A{u2))A{u3)) 

J(a:,oo)" 

...{{{I + R,)A{Un))A{ui)) 

= {-rr{{I + RM{u)Mu)Y. 

□ 

The identities in the following Lemma can be found in or deduced from Tracy- 
Widom Ii26il. 



Lemma 9.3. 

(9.1) R,{x,x) = {{I + R,)A,A) 

{9.2)(^ + -^ + -f]R,{u,v) = -((/ + /?,)A(m))((/+/?,)A(v)) 
\ax du dv J 

(93)(-^ + -^\l+R,)A{u) = {l + R,)A'{u)-{I + R,)A{u){{I + RM,A). 
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(9.4) 2((/ + /?,)A',A) - {{I + R,)A,Af = - {{I + R,)A{x)f 

(9.5) 2{{I + R,)A,A") - {{I+R.M',A') = x{{I+R,)A,A) 

(9.6) ^{I+R,)A'iu) = -R,{u,x){I+R,)A'{x) 



(9.7) —{I + R,)A'{u)=u{I + R^)A{u)-2{{I + R,)A',A){I + R,)A{u) 
au 

+ {{I+R,)A,A){I+R,)A'{u)+R{u,x)iI+R,)A'{x). 
Lemma 9.4. We have 

Q'oix) = {{I + R,)A,A) and ^(.{x) = -{{I + R,)A{x)f 



1 



(2[,"(x) = -((/+/?,)A(x))((/+/?,)A'(x)) + ((/ + /?,)A(x))2((/+/?,)A,A). 



Proof. One computes 



d 

dx 
d 

dx 



/poo 1 pOO /»oo \ 

/ K^^\u,u)du + -j K''^^\u,w)K'^^\w,u)dudw ^ 

Kq{u,w\)Kq{w\_,W2)K^^'' {w2,u)dudwidw2 + ■■■ 



noa poo poQ 



\+3 

/•CO 

K^^\x,x)+ I K^'^\x,v)K^^\v,x)dv 

X 

is:*") {x, wi (wi , W2)is:(°) (w2,x)^i^wl + . . . 



+ 



= R,{x,x) = {{I + R,)A,A), using ^ 

by the Neumann series for Rx = {I — K^^^)^^ — I. Moreover, using the previous 
result, 

///s d , , ( d d d \ 

e'o'(.) = -/e.(x,.) = ^ + ^ + ^U.(«,v) 



-{{I + R,)Aiu))iI + Rx)A{v)) 
-{{I + Rx)A{x)y-, 



using (9.2 1 
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and again using the result just obtained and using identity (|9.3[), 



= -2{{I + R,)A{x))-{I+R,)A{x) 

= -2{{I + R,)A{x)) + (/ + /?.)A(«) 



= -2((/ + /?,)A(x))(/ + /?,)A'(^) + 2((/+/?,)AW)'((/ + /?.)A,A), 
proving Lemma [94] 
Lemma 9.5. 

21 (x) =r((/ + /?,)A,A) 

Qiix) = -^((/ + /?,)A(x))2 
^3 

eaW = -{{{I + R,)A{x)f{{I + R,)A,A)-{I + R,)A{x){I + RM'{x)) 



□ 



rx 



((/+/?,)A,A). 



Proof. Indeed, by (|7^ and (|2J3), 
(9.8) !2i(x) = -TrLi = 

= r/ du A{u) [ A{u) + I Rx{u,v)A{v)dv 



- Tr kI'^ -Tr R^K^'^ 



poo 

r / A{u){{I + R)A){u)du = r{{I + R^)A,A) 

Jx 



Computing Q2 by (7.8 1 involves TrL2 and TrLj. Since K2\u,v) has a symmetric 
and skew-symmetric part, and since I + Rx is symmetric, and remembering the 
form of K2\u,v) in (2.13 1, we have (by symmetry) and the fact that a symmetric 
times a skew-symmetric operator is traceless. 



TtL2 = Tt{I + Rx)K: 



--Tr(/ + /?.,)(A'(M)A(v)+A(M)A'(v)) 
-r^{{I + Rx)A',A). 



Hence, combining this result with the computation of TiLj in Lemma 
finds, using (|9.4[), 



9.2 



one 



(9.9) !22 = -Tr(L2 + 



-{2{{I + Rx)A',A)-{{I + Rx)A,A) 



-((/ + /?,)A(x))2. 
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The computation of ^3, by (7.8 1, involves TrLs, TrLiL2 and TrLj. Using again the 



fact that a symmetric times a sk ew-symmetric operator is traceless, one reads off 



(r) 

from the form of (see (|2.13|)) the following: 



TrL3 = Tr{I + R,)K';' = --{{{I +R,)A,A") + {{I +R,)A' ,A')) 



3 



/^{2{{I + R,)A,A")-{{I+RM',A')) 



and using 

{I+R,)K^^ = -^(((7 + /?,)A'(«))A(v) + ((7 + /?,)A(m))A'(v)) 

-L[{{i+r;)A'{u))a{v)-{{i+rM{u))a'{v)). 



one computes 



TrLiL2 = 1:v{1 + R^)k'(\i + R,)K^^ 

.-3 



^ // {{I+R,)A){wi)A{w2m + R,)^){w2)A{w,)dw,dw2 

+ "^11 {{I + RM){Wi)A{w2){{I +Rx)A){w2)A'{wi)dWidW2 

J J {{I + R,,)A){wi)A{w2){{I + RM'){w2)A{wi)dwidw2 



r2 



2 

' 2 

3 / {{I + R,)A,A){{I + R,)A',A) 



J I ((/ + /?,)A)(W1)A(W2)((/ + /?,)A)(W2)A'(W1)JW1JW2 



+{{I + R,)A,A){A',{I + RM) 
2 ( {{I+R,)A,A){{I + R,)A',A) 

-{{I + R,)A,A){A',{I+R,)A) 

r3((/ + /?,)A,A)((/+/?,)A',A). 



r 
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(9.10) 23 = -Tr(L3+LiL2 + ^L3 



Putting the pieces together and using (7.8 1 and using Lemmas 9.2 and 9.4 one 
obtains 

1 

— i 
3 

^3 / ((/ + /?,)A,A") + ((/ + /?.r)A',A') 
-3((/+/?,)A,A)((/+/?,)A',A) 
+((/ + /?,)A,A)3 

+L{2{{I+RMA")-{{I+RM'A')) 

3 / {(l + R^)A,A") + {{I + RM'A') 
L -{(l + R^)A,A){{I + R,)A',A) 

^ y +{{I + R,)A,A)(^{{I + R,)A,Af-2{{I + R,)A',A] 
+ ^(2((/ + /?,)A,A") - {{I + R,)A',A')) 



-{I+R,)A{x){I + R,)A'{x) + {{I + R)A{x))^{I+R,)A,A)) 



rx 



+ -{{I + R,)A,A), 



using in the last equality ( |9.4| ), ( |9.5| ) combined with ( |9.1 1| ) below. Then, using the 
differential equation uA{u) = A"(m), one checks: 

(9.11) 

((/ + /?,)A',A') 

= -A(;c)(/ + /?;i-)A'(x) - i^il + RM' ,A) (by integration by parts) 

au 

= -A{x){I + R^)A\x) - {{I+R^)A{u),uA{u)) + {{l + RM\A){{i + Rx)A,A) 

-{R{u,x),A{u)){I + R)A'{x), using 
= -(/+/?,)A(x)(/ + /?,)A'(x)-((/ + /?,)A,A") + ((/ + /?,)A',A)((/ + /?,)A,A). 

□ 



Proof of Proposition \9. 1 \ The formulae follow immediately from comparing the 

□ 



formulae of Lemmas 19.41 and B. 
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